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The pair correlation function of an expanding gas is investigated with an emphasis on the BEC- 
BCS crossover of a superfluid Fermi gas at zero temperature. At unitarity quantum Monte Carlo 
simulations reveal the occurrence of a sizable bunching effect due to interactions in the spin up- 
down channel which, at short distances, is larger than that exhibited by thermal bosons in the 
Hanbury-Brown and Twiss effect. We propose a local equilibrium ansatz for the pair correlation 
function which we predict will remain isotropic during the expansion even if the trapping potential 
is anisotropic, contrary to what happens for non-interacting gases. This behavior is understood to 
be a consequence of the violation of scaling of the pair correlation function due to interactions. 
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Recent experimental studies of two-body correlations 
in an expanding atomic cloud are opening new perspec- 
tives in the study of quantum-statistical effects in ul- 
tracold gases. These include the measurement of the 
Hanbury-Brown and Twiss effect [J] in a dilute Bosegas, 
both concerning spatial and temporal correlations [^ |3| , 
measurements in the Mott insulating phase of ultracold 
gases in an optical lattice [J|, as well as the study of 
atomic pair correlations in a Fermi gas after dissocia- 
tion of molecules near a Feshbach resonance Q ■ In some 
cases the measured quantities are related to the real space 
properties {i.e. to the pair correlation function), whereas 
in other cases they refer to correlations in momentum 
space. In all cases, measurements are done after a period 
of free expansion since the gas in situ does not provide 
sufficient optical resolution. 

The pair correlation function is defined as 



(2). s n(2) ri,cr;r2,cr' 

5^j'r 1,1-2 = —, ^-, jr , 

n(ri,cr)n(r2,cr') 



(1) 



where 



i(^^(ri,(T;r2,cr') 



{i^t{n)^lir2)i^.,{r,)Mri)) , (2) 



is the second order correlation function and n(r, a) is the 
density of the a spin component. The spin indices cr, 
a' are used if there is more than one species of atoms 
present. The pair correlation function is a key quantity 
in many-body physics, being sensitive to both statisti- 
cal and interaction effects. The statistical effects appear, 
for example, in non-interacting gases of identical bosons 
above the critical temperature, where they are respon- 
sible for a characteristic bunching effect which has been 
recently the object of experimental measurements |2]. 

In this Letter we discuss the case of a Fermi gas, with 
equal population of two spin states, close to a Feshbach 
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resonance where the value of the scattering length a can 
be tuned by varying the external magnetic field. The op- 

(2) 

posite spin gl. correlation function is strongly affected 
by interactions and for a homogeneous system has been 
calculated using fixed-node diffusion Monte Carlo sim- 
ulations 6]. The result at T = is shown in Fig. Q] 

where g\/ (s) is plotted as a function of the dimensionless 

variable kps for various configurations. In the inset we 

plot the integrated quantity N-^ils) — n f^ iirr^drgiy (r) 

giving the average number of atoms of opposite spin 

within a sphere of radius s around a given atom. Here 

n — kp/{6TT'^) is the single spin density and kp is the 

Fermi momentum. On the BEC side of the resonance 

{kpa = 0.25, dashed line) there is a pronounced bunching 

effect at short distances due to the presence of diatomic 

molecules and the pair correlation function rapidly ap- 

r2i 
proaches the uncorrelated value 5} i (s) — 1 (thin solid 

line) at distances of the order of the scattering length. As 
a consequence, the integrated probability N^ils) reaches 
1 at kps ~ 0.25 in contrast to the uncorrelated gas where 
the same value is reached much later, at fc^^s = 2.4. In 
the unitary limit {1/kpa = 0, solid line) g\. is a univer- 
sal function of kps, the inverse of kp providing the only 
length scale of the problem. In this case the bunching ef- 
fect is reduced as compared with the BEC case, but is still 
important (for example at kps = 1 the value oi N'^i{s) is 
0.5 whereas the non-interacting value is negligible). For 
comparison we plot the pair correlation function for a gas 
of non- interacting bosons with density n — fc|./(67r^) at 
the Bose-Einstein critical temperature Tc (dotted line). 
At short distances g^^\s) = 2 and the corresponding in- 
tegrated probability N{s), while larger than that of the 
uncorrelated gas, is still quite small compared to the uni- 
tary case [3. The sizable effect exhibited by the g\\ (s) 
correlation function in interacting Fermi gases is due to 
the presence of correlations resulting in a 1/s^ behavior 
as s —5- [3 which characterizes not only the deep BEC 
limit, but the whole crossover Q . The short range behav- 




FIG. 1: (color online) . Spin up-down pair correlation function 
g[,{s) of a homogeneous system. Dashed line (red online): 
Fermi gas in the deep BEG regime {l/kpa = 4); solid line 
(blue online): Fermi gas at unitarity {l/kpa = 0); dotted line 
(green online): ideal Bose gas with density n — fc|-/(67r^) at 
the Bose-Einstein critical temperature; thin solid line (black 
online): uncorrelated gas with g (r) = 1. Inset: Integrated 
pair correlation function A'^i4.(s) for the same configurations. 
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lor of (^li has recently been the subject of experimental 
studies using spectroscopic techniques [lOj. At unitarity, 
the s —f limit of s'^gi, (s) is fixed by many-body effects 
and the quantum Monte Carlo simulation yields the value 
(kpsyg^Ms) -^ 2.7. In the BEC limit one finds instead 

the value (fcFs)^.9t| (s) -^ 37r/{kpa) characterizing a free 
molecule. 

As already mentioned it is difficult to directly mea- 
sure the pair correlation function in real space while 
the atomic gas is confined in the trap since the cloud 
is too small to optically resolve features at interatomic 
distances |Tl|. The gas needs to be released, undergoing 
an expansion, before any measurement can be made. If 
the scattering length is kept constant during the expan- 
sion, the calculation of the time dependent evolution of 
the pair correlation function is strongly affected by in- 
teractions and represents a challenging theoretical prob- 
lem for which few exact solutions are known. One of 
them is the unitary Fermi gas released from an isotropic 
harmonic potential ^12, 13] . In the more general case 
(including anisotropic traps or systems which are not at 
unitarity) the behavior of the expansion can be neverthe- 
less predicted if the system is in the hydrodynamic (HD) 
regime, corresponding to a local equilibrium condition. 
In this case one can use the ansatz 



(2) 



g^'>{ri,r2,t)=gi:jjs,nir,t)), 



(3) 



for the pair correlation function, based on a local density 
approximation. Here r = (ri -I- r2)/2 and s = ri — r2 are. 



respectively, the center of mass and relative coordinates 

(2) 

and 5jjom ^® ^^^ P^^^ correlation function for a homoge- 
neous gas with density n(r, f). From an experimental 
point of view the ansatz © relates the locally measured 
density n(r, t) to the locally measured g^'^\ From a theo- 
retical perspective it gives us the value of g^^'' during the 
expansion if we know the equilibrium form of gl.omi^) 
(from, e.g., quantum Monte Carlo calculations) and the 
time evolution of the local density n(r, t). In particular, 

at unitarity, where g/^J^ depends only on the combina- 
tion kps, the expansion acts like a microscope, the value 
of kp ex n^'^ being reduced as a function of time. 

The dynamics of the density n{r,t) can be analysed 
employing the HD theory. This theory has been quite 
successful in predicting the behavior of n(r,i) in sev- 
eral important configurations. For example in the case of 
a dilute Bose gas the hydrodynamic predictions for the 
time evolution of the density agree in an excellent way 
with the exact solutions of the Gross-Pitaevskii equation 
and with experiments (see, for example, Ref. [Ijl). They 
have been used to study, among others, the collective os- 
cillations 1151 a. nd the expansion after release from the 
trap ^^ [13, [i3 ■ If the trapping has the harmonic form 
^(i") = J2i=x y z '^'^?^?/2 and the chemical potential fol- 
lows a power law dependence on the density, fJ-(n) ex n'"' 
the HD predictions for the expanding density profile can 
be written in the scaled form 



n{r,t) 



1 



bxbybz 



no{x/bx,y/by,z/bz) 



(4) 



with the scaling coefficients 6,; obeying the system |lq 



6,:- 



b^{bxbyb,)^ 



0, 



(5) 



of nonlinear coupled equations with the initial conditions 
bi{t = 0) = 1 and no(r) = n{r,t — 0). The polytropic 
dependence of the chemical potential applies both to the 
case of a dilute Bose gas (7 = 1) and to the unitary 
limit of a Fermi gas (7 = 2/3). When released from 
an anisotropic trap, the aspect ratio of the expanding 
density profile exhibits an inversion from the cigar to 
the disc geometry (and viceversa) as a consequence of 
the HD forces which are more active along the directions 
where the gradient of the density is larger. Therefore the 
density profile at long times will be anisotropic. On the 
other hand we see from Eq. Q that the pair correlation 
function will keep its isotropy as a function of s. 

The ansatz 0, which assumes local equilibrium for 
the description of the expansion, requires more strin- 
gent conditions compared to the usual local density ap- 
proximation applied to stationary configurations. In the 
latter case, the only relevant condition is that the typi- 
cal length scale, where g^'^^ approaches the uncorrelated 
value g'^' = 1, be much smaller than the size of the cloud. 
The use of Eq. (Q , instead, also requires that adiabatic- 
ity be ensured during the expansion. In the superfluid 
regime, HD theory provides a justification for the local 



equilibrium approximation. We might however wonder 
whether superfluidity is preserved during the expansion. 
In a Fermi gas at T = 0, the disappearance of superfluid- 
ity could take place through pair breaking mechanisms. 
Since a large number of quantized vortices are observed 
after expansion, recent experiments [l^ suggest that the 
system remains superfluid during the expansion on the 
BEC side of the crossover. The robustness of superflu- 
idity is also suggested by the fact that the gap is much 
larger than the typical oscillator frequency (whose inverse 
is the timescale of the expansion) and that, as a conse- 
quence, pairs cannot break during the expansion. This 
argument holds also at unitarity where the gap is of the 
order of the Fermi energy. On the other hand, super- 
fluidity is more fragile on the BCS side because the gap 
becomes exponentially small during the expansion (for a 
measurement of the pairing gap in ultracold Fermi gases 
along the BEC-BCS crossover see Ref. pw|). 

The experimental measurement of the pair correlation 
function and the verification of Eq. (O requires that 
certain conditions be met. In particular, the expansion 
time should be long enough for features at the inter- 
atomic level to be resolved. By choosing, for example, a 
trapped gas with central density equal to n ~ 10^'^ cm~^ 
(corresponding to an inverse Fermi momentum equal to 
kp ~ 100 nm) and isotropic trapping, we find that, after 
a time equal to wi ~ 20, the density is reduced to a value 



10 cm '' corresponding to k^ 



3/im. This value 



is compatible with the present limits of optical resolu- 
tion 21]. At the same time, if we choose a large value of 
the scattering length the system will keep the condition 
of unitarity kpa ^ 1 also at the observation time. 

The HD description of the expansion points out the 
occurrence of important qualitative differences in the ex- 
pansion of interacting and non-interacting gases (see Ta- 
ble P). Contrary to what happens in the HD regime, 
when released from an anisotropic trap the density of a 
non-interacting gas (e.g., fermions or uncondensed ther- 
mal bosons) becomes isotropic at long times during the 
expansion, reflecting the isotropy of the initial momen- 
tum distribution J22l | . The pair correlation function also 
exhibits drastic differences. In fact, in a non-interacting 
gas the dependence of g^^^ on the relative coordinate s 
will undergo an anisotropic expansion. For example, for 
non-interacting bosonic or fermionic gases trapped in a 
harmonic potential at temperatures much higher than 
that of quantum degeneracy, the pair correlation func- 
tion can be shown to exhibit the following behavior |23| 



J2) 



(s,i) = 1 ± exp 



-27r 



E 



\Ucjft^ + l) 



(6) 



where At — \/2TTh?/mkBT is the thermal de Broglie 
wavelength. Here the plus sign is for bosons, revealing a 
bunching effect (as discussed previously), while the mi- 
nus sign is for fermions which exhibit anti-bunching as 
a consequence of the Pauli principle. For large times. 



TABLE I: Role of anisotropy in the expansion of the density 
and of the pair correlation function of a trapped gas. 



Behavior after expansion 


n{r,t) 


ff^^'(s,t) 


Ideal gas 

(bosons T > Tc or fermions) 


isotropic 


anisotropic 


Hydrodynamic regime 


anisotropic 


isotropic 



the typical value of the decay length of these statistical 
correlations depends on direction and scales as Xrujit = 



y/ATrht/{mRi), a value that has recently been used to 
analyze the data in the experiment of J3| with thermal 
bosons. Here Ri = ^2kBT jraiJl is the in situ radius of 
the thermal cloud in the z-th direction. 

Let us finally discuss the problem of the expansion 
from the general point of view of scaling, characterized 
by a coordinate transformation of the type: 



X -^ x/b.^{t), y ~* y/by{t), z -^ z/b^(t) 



(7) 



An exact scaling solution is known to be obeyed by the 
many-body wave function of a unitary Fermi gas released 
from an isotropic harmonic potential |l2l | . In this case all 
correlation functions exhibit the same scaling behavior 
and one finds ri(r,t) = na{r/b)/b^ and n'^^^(ri,r2, i) — 
n(2)(ri/6,r2/6,0)/66, with b = VUWTT. The HD pic- 
ture reproduces this behavior, in fact Eqs. Q)-®, with 
uji = Lu, bi(t) = b{t) and 7 — 2/3, yield the same time de- 
pendence for the density profile. Agreement is also found 
for the time dependence of the pair correlation function. 
In fact at unitarity the correlation function evaluated for 
a homogeneous medium depends only on the combina- 
tion kps oc 'n}/'^s so that Eq. Q can be written in the 
scaled form 



9 






ny\v/b)slb 



(8) 



which, for isotropic harmonic trapping, agrees with the 
scaling result ^(^^(ri, r2, t) = g'''^\vi/b,V2/b,Q). Nat- 
urally the agreement holds only for values of s much 
smaller than the size of the system where the local den- 
sity approximation can be properly applied. 

The applicability of scaling to describe the expansion 
of the pair correlation function is however not general. In 
fact for an interacting gas the boundary conditions on the 
many-body wave function impose restrictions which are 
in general incompatible with scaling. At distances much 
shorter than the interparticle separation, but much larger 
than the effective range of interatomic interactions, the 
wave function is determined up to a proportionality con- 
stant by the Bethe-Peierls boundary condition: for any 
value of the scattering length a and for all pairs of atoms 
ij, the many-body wave function ^ obeys the condition: 



*(|r. 



0)cx 



1 
a 



(9) 



where the hmit is taken fixing the position of all other 
atoms as well as that of the center of mass (r^ + Yj)/2. 
This boundary condition has two characteristic features: 

I) it is isotropic and II) it introduces a length scale a. 
Any scaling solution of the form of Eq. ((TJ will run into 
two kinds of difficulties. I) requires that all the hi{t) be 
identical, which is in general incompatible with the free 
expansion from an anisotropic trap. II) means that, even 
if the trap is isotropic, leading to identical bi{t), a rescal- 
ing of the coordinates would change the ratio a/\vi — r^l 
which is required to be fixed. These problems do not ap- 
pear when a = or a ^ oo, i.e., in the non- interacting 
and isotropic unitary gas |2J]. The above argument ap- 
plies to 5^^-' and to all higher order correlation functions 
as they all obey the boundary condition ©. However, it 
does not affect the density n{r,t) for which scaling laws 
do exist, e.g., for HD systems. 

From I) we can now better understand the behaviour 
of f/'^^-' summarised in the Table. It implies that g'^^ 
must remain isotropic (at least for small |ri — r2|) for 
all interacting gases whereas in the non-interacting limit 
this constraint is removed and g^'^' is allowed to exhibit 
the anisotropic scaling of Eq. ©. On the other hand, 

II) leads to a violation of scaling which dramatically il- 
lustrates the difference between the BEG and unitary 
regimes: the latter has perfect scaling of gl/ if the gas 



is isotropic. Conversely, in the molecular gas in the deep 
BEC regime, although the density n(r, t) scales during 
the expansion as in ordinary Bose-Einstein condensates, 

(2) 

5+ 1 does not since it now contains information about the 
molecular bound state whose size is fixed by the scatter- 
ing length a and is independent of the local value of the 
density. 

In conclusion, we investigated the spin up-down pair 
correlation function of a Fermi gas at T = in the BEC- 
BCS crossover. Under the condition that local equilib- 
rium is ensured during the expansion, we addressed the 
non trivial problem of the time dependence of the corre- 
lation function after release of the gas from the confining 
potential. In contrast to the behavior of non-interacting 
gases, the pair correlation function violates the scaling 
solution (O and is predicted to remain locally isotropic 
during the expansion even for anisotropic trapping. At 
unitarity, the absence of other relevant length scales be- 
side kp makes the expansion act like a magnification 
lens as the average distance between particles becomes 
visible. 
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